Reflections on dissipation associated with thermal convection 
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Buoyancy-driven convection is modelled using the Navier-Stokes and entropy equations. It is 
first shown that the coefficient of heat capacity at constant pressure, Cp, must in general depend 
explicitly on pressure {i.e. is not a function of temperature alone) in order to resolve a dissipation 
inconsistency. It is shown that energy dissipation in a statistically steady state is the time-averaged 
volume integral of — and not that of — aT-^. Secondly, in the framework of the anelastic 
equations derived with respect to the adiabatic reference state, we obtain a condition when the 
anelastic liquid approximation can be made, 7 — 1 << 1, independent of the dissipation number. 

PACS numbers: 44.25. -|-f, 47.10. ab, 91.45.Ga, 91.45.Kn 



Convection in the geophysical context is often associ- 
ated with compressible effects (as early as in the work of 
Carnot [l[) and simultaneously with a significant energy 
dissipation , of the same order of magnitude as the heat 
flux. In the case of convection in the Earth's mantle, it is 
important to determine the distribution of dissipation as 
it is strongly coupled with effective viscosity and plays a 
role in the structure of convection Q , while the amount 
of dissipation in the liquid Earth's core is a direct mea- 
sure of the energy available to the geodynamo f2] . In this 
paper we shall be interested specifically in dissipation in 
'liquids', defined as those fiuids with a small product aT 
compared to unity, where a is the coefficient of thermal 
expansion at constant pressure and T is the thermody- 
namic temperature. This product is unity for perfect and 
semi-perfect gases, but is usually very small for liquids, 
of order 0.05 both in the outer core of the Earth, and in 
the mantle which are hence both usually considered as 
'liquids' from the perspective of convection. 

Let us consider a simple configuration of a region 
of Newtonian fluid with no internal heat production, 
bounded with solid walls or stress-free boundaries and 
prescribed temperatures or heat fluxes within a uniform, 
constant gravity field g. It might be a simple Rayleigh- 
Benard configuration with an imposed temperature dif- 
ference. It might also be more complex, we simply require 
that the boundary exerts no work on the fiuid region and 
that a statistically steady state can be reached so that 
the net heat input is zero. Density p, velocity u, pressure 
P, temperature T and entropy s satisfy the continuity, 
Navier-Stokes and entropy (or heat transfer) equations 
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where (p — —kSJT is the conduction heat flux (with ther- 
mal conductivity k), r is the stress tensor and e is the 
tensor of the rate of deformation. Rewriting ([3]) using 



the thermodynamic relation Tds = CpdT — (well- 
known equation (jA.37l) re-derived in the appendix), with 
Cp the specific heat capacity at constant pressure, leads 
to 
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A classical expression for the energy dissipation is ob- 
tained by integration of equation @ over the whole fluid 
domain, under the tentative assumption that Cp can be 
taken as uniform, and by taking its time-average (^Q^ 
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where (.) denotes time-averaged volume integral. How- 
ever, the integration of the dot product of the Navier- 
stokes equation ([2]) with velocity u, assuming no work is 
done through the boundary, leads to 



(e : t) = - (u • VP) = 



(6) 



which is very different from equation ([5]) for liquids as 
aT « 1. The inconsistency is resolved when it is rec- 
ognized that Cp can neither be uniform nor a function 



of temperature only. As Cp 
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by definition, with 



H the enthalpy, the term pcp-^ in equation ^ can be 
integrated as follows 
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The first term of the right-hand side is always zero (this 
would be true for any function of state, not just for iJ) 
in the statistically steady case. The second term is ex- 
pressed using equation (IA.38|) of the appendix, so that 
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By making the hypothesis of a uniform Cp, an important 
term has been forgotten, thus invalidating expression ([5]). 
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The correct averaging of the heat equation, just hke the 
balance of mechanical energy, leads to (|6]). Incidently, 
we may just write two other exact expressions for dissi- 
pation, at least as fundamental as (jS]), obtained from the 
integration of (U) where entropy variations are not ex- 
panded in terms of pressure and temperature variations: 
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Equation ([5]) will now be used, within the anelastic 
approximation (AA) Q, to determine a condition when 
the anelastic liquid approximation (ALA) [1, 0] can be 
made. We assume that convection is sufficiently vigor- 
ous to drive the fluid close to an iso-entropy state, com- 
monly called 'adiabatic'. All quantities will be expanded 
around their adiabatic reference state denoted with the 
subscript while fluctuations are indicated by primes. 
The reference state is dependent on height only and is 
defined such that the entropy sq is uniform and that the 
hydrostatic equation is satisfied dP^/dz = —pog, where 
z is the vertical upward coordinate. Equation (jA.37[) 
shows that, when entropy is vigorously mixed, a relation 
develops between pressure and temperature variations. 
In addition, when hydrostatics is considered, this leads 
to the so-called 'adiabatic gradient' 
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The main purpose of the anelastic approximation is to 
eliminate sound-waves by replacing the general continu- 
ity equation ([T]), by its zeroth-order expansion 

V • (pou) = 0. (11) 

The first order expansion of the entropy equation (|3]) is 
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Writing the first term in a conservative form introduces 
a advection term for the adiabatic temperature profile 
which is expressed using the adiabatic gradient (first term 
on the right-hand side of the following equation) 
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The derivation of the anelastic momentum equation is 
then based on the expansion of density fluctuations in 
terms of entropy and pressure fluctuations 
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Using the Maxwell relation associated with ()A.34[) and 
the expression for the adiabatic temperature gradient ob- 
tained from (jA.37l) for the first term, and hydrostatics for 
the second, one obtains 
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Substituting in the Navier-Stokes equation ([2]) and se- 
lecting the lowest order terms leads to 

Po-p^ = -yP^ : s ez + — -^P Bz+y-T, (16) 
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where is the vertical unit vector. The first and third 
terms on the right-hand side can be put together in a 
single conservative term 
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From this point, we can derive important consequences 
in terms of energy dissipation. In (fTTl) the work of 
— poV (P'/po) integrated over the volume is zero. There- 
fore the mechanical works of the first (—VP') and third 
(P'Vpo/po) terms of ([TCI) are opposite to each other. 
As the expression of the dissipation — (DP/Dt) becomes 
- (DP'/Dt) = - (u • VP'), we have 
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where Kg is the incompressibility at constant entropy. 
Introducing c = y^Ks/ p, the celerity of sound waves, one 
can finally express the energy dissipation in terms of the 
correlation of velocity and pressure fluctuations 
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However, energy dissipation must also be equal to the 
work done by the second term on the right-hand side 
of (|16|) . in agreement with the expression of dissipation 
associated with ([T3l) and the general expression 
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now in terms of the correlation between velocity and en- 
tropy fluctuations. Expanding entropy fluctuations in 
terms of temperature and pressure fluctuations, using 
(|X37)) . leads to 



(e : r) = {aopogu.T ) - ( u,P ) 

\ CpO / 
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Comparing equations (IT^ and (|21l) shows the relative 
importance of the pressure and the temperature terms in 
d^I]) . The coefficient within the pressure/ velocity corre- 
lation term can be written 
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where 7 = Cp/cy is the ratio of specific heat capacities 
at constant pressure and at constant volume. The sec- 
ond equality in (j22p can be retrieved from the general 



3 



Mayer relation Cp - = -T/ p'^{dP/dT)\p{dp/dT)\p . 
Under the assumption that velocity and pressure fluctu- 
ations are reasonably correlated (i. e. the time- averaged 
integral of u^P' on horizontal surfaces is positive at all 
heights, or at least that this condition is not too much 
violated), the contribution of pressure/velocity correla- 
tions to dissipation in (1211) can be bounded as follows 



following relation 



/mm {e:T)< 



u,P' ) < f^^ax (e-.T), (23) 



where fmin and fmax denote the minimum and maxi- 
mum respectively of the quantity 70 — 1 over the height 
of the fluid domain. Of course, in a special case where 
that quantity would be uniform, we would have the ex- 
act relationship {^{a^Tag / Cpoju^P') = (70 — 1) (e : r). It 
is hence safe from the point of view of energy dissipa- 
tion to express entropy fluctuations s' (see I20p in terms 
of temperature fluctuations only T' (see I23p in equation 
(fT7)) under the condition 



1-1 = << 1, 



(24) 



constituting the proper condition of application of the 
anelastic liquid approximation. While entropy fluctua- 
tions are replaced by temperature fluctuations in ([T7| . 
CpoT' must be substituted to Tqs' in (|13p to ensure a con- 
sistent energy balance. When p4|) is satisfied, the follow- 
ing anelastic liquid approximation equations are derived 
from the anelastic approximation (fTTj) . (fT3|) and (flT)): 
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Let us mention here, for completeness, that a so-called 
truncated anelastic liquid approximation (TALA) has 
been used (for instance in (lOi] ) in which the contribu- 
tion of pressure fluctuations to density fluctuations are 
also neglected in equation (ITSI) . As a consequence, the 
term —poV {P' /po) in equation (|27)) above is changed for 
—VP'. However, as pointed out in [11], this change in- 
troduces an imbalance between energy dissipation calcu- 
lated from the dynamical equation and heat dissipation in 
the thermal equation. Such a TALA formulation should 
be avoided when the calculation of energy dissipation is 
an issue. 

Our condition (1241) is in contrast to a scaling law usu- 
ally derived [1, Il2| from the anelastic liquid approxi- 
mation momentum equation p7p expressing the typical 
magnitude of pressure fluctuations as a function of tem- 
perature fluctuations, P' ^ a^pogHT' , where H is the 
typical length-scale of the fluid domain. This leads to the 
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which would imply that the anelastic liquid approxima- 
tion is valid when (ao7o)I? << 1, where the dissipation 
number is _D = aogH/cpo. This heuristic scaling law ([25)) 
involving D, hence H and g, is not compatible with the 
exact condition ([M]) that is expressed solely in terms of 
fluid properties. Actually, from equations ([TO)) and ([^ . 
one can infer a more adequate scaling of the pressure 
fluctuations than P' aopogHT' , which is valid in the 
general anelastic approximation (liquid or not) 



P' - KtoaoT', 



(29) 



where Kt = K^/"/ is the incompressibility at constant 
pressure. 

Let us now discuss the implications of the results de- 
rived in this paper. Concerning the 'exact' model ([1]), 1^ 
and ([3]) , we have shown that it is crucial to use thermody- 
namically consistent thcrmo-physical coefficients. Specif- 
ically, considering a fluid with constant heat capacity Cp 
and small product aT of temperature by thermal ex- 
pansion coefficient leads to a severe inconsistency when 
determining the energetic dissipation. Whatever aT, dis- 
sipation is always equal to {DP/ Dt) . Once compressibil- 
ity effects have been taken into account in an anelastic 
model, such an inconsistency is no longer to be feared 
and one can safely use constant uniform values for Cpo 
and other coefficients. Using the correct expression for 
dissipation, we have been able to determine rigorously 
the condition 7 — 1 < < 1 for the anelastic liquid approx- 
imation to be valid. That condition is only dependent 
on materials properties, not on gravity or length-scale. 
Note that the condition ([^ can be expressed with the 
Griineisen parameter, F = aKs/{pCv) used in Solid state 
Physics, as (?o?T jcp = 7 — 1 = r7aT. The Griineisen 
parameter is close to unity in the mantle and in the 
Earth's core and so is the ratio of heat capacities. This 
is the reason why the condition ([M)) is close to aT << 1 
in practice concerning the dynamics of the deep Earth. 
For the Earth's liquid outer core and for the mantle, the 
dimensionless number in (|24|) is of order 0.05 justifying 
the anelastic liquid approximation. 

From the correct expression for dissipation ((51), we have 
obtained equations ([^ and ([?7|) . in the anelastic liquid 
approximation, which both lead to the following expres- 
sion for dissipation 



{aopogT'uy_ 



(30) 



Surprisingly, Hewitt et al. [5| reach the same conclusion 
when starting from their expression ([5]): 
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This expression is indeed approximately correct for liq- 
uids and the apparently very different expressions ([6]) and 
([5]) are nearly equal in the anelastic liquid approximation 
(DP/Dt) « (aTDP/Dt). This is due to the large cancel- 
lations in ^ which allow the same result to be reached 
when multiplied by a small term aT, with appropriate 
fluctuations (here aoT' to the leading order). We only 
expect significantly different values from expressions ^ 
and ^ when aT is neither too small nor too close to 
unity, where ([6]) and ([5]) become similar. 

To conclude, we hope that we have contributed to a 
clarification of the expression of dissipation in a con- 
vective system and of the anelastic liquid approxima- 
tion. We have not systematically investigated the conse- 
quences of these findings in the different possible fields 
of application: dynamics of the Earth's core and man- 
tle, dynamics of giant planets, ice planets, super-Earth 
exoplanets. Regarding the application of the anelastic 
liquid approximation for the terrestrial planets obeying 
a generic Murnaghan equation of state |13|, we have no 
indication that the new criterion ([M)) . 7 — 1 << 1, dif- 
fers significantly from the classical criterion aTD << 1. 
We have considered the case of silicate planets from one 
up to possibly ten Earth masses. The parameter 7 — 1 
is always less than 0.08, and that value decreases very 
quickly away from the surface to much lower values in 
the bulk of the mantle (decreasing also when planetary 
mass increases), so that the anelastic liquid approxima- 
tion is indeed well justified for the mantle dynamics of 
the Earth and even better justified for super-Earths. 



Useful thermodynamic relations 

From the Gibbs equation expressed in terms of the 
Gibbs free energy G 
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the following Maxwell relation is derived, expressing the 
partial derivative of entropy with respect to pressure at 
constant temperature 
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The partial derivative of entropy with respect to temper- 
ature at constant pressure is obtained as follows. Two 
expression for dH are first written, one from Gibbs rela- 
tion, the other from the definition of Cp 

1 



dH = TdS+-dP, 
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The difference of these expressions above leads to 
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This shows that the partial derivative of entropy with re- 
spect to temperature at constant pressure is Cp/T. Along 
with equation (|A.33|) . this leads to 
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Comparing equations (jA.36|) and (jA.37P leads to an im- 
portant expression used in this paper 
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The partial derivative of this last equation with respect 
to temperature at constant pressure leads to expressions 
of the derivative of Cp with respect to pressure at constant 
temperature 
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